Evidence of short time dynamical correlations in simple liquids. 
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We report a molecular dynamics (MD) study of the collective dynamics of a simple monatomic 
liquid -interacting through a two body potential that mimics that of lithium- across the liquid-glass 
transition. In the glassy phase we find evidences of a fast relaxation process similar to that recently 
found in Lennard- Jones glasses. The origin of this process is ascribed to the topological disorder, 
i.e. to the dephasing of the different momentum Q Fourier components of the actual normal modes 
of vibration of the disordered structure. More important, we find that the fast relaxation persists in 
the liquid phase with almost no temperature dependence of its characteristic parameters (strength 
and relaxation time). We conclude, therefore, that in the liquid phase well above the melting point, 
at variance with the usual assumption of un-correlated binary collisions, the short time particles 
motion is strongly correlated and can be described via a normal mode expansion of the atomic 
dynamics. 



PACS numbers: 63.50.-Hx, 64.70.Pf, 61.43. Fs, 62.60.-hv 



I. INTRODUCTION 

Many experimental investigations by Inelastic Neutron 
(INS) and X-ray (IXS) Scattering techniques, as well as 
different theoretical and numerical studies have been de- 
voted since more than twenty years to the understanding 
of the collective dynamics of simple liquids. In partic- 
ular, the efforts have been focused on those dynamical 
ranges that go beyond the two limiting cases of single 
particle (high exchanged momentum, Q) and hydrody- 
namic (low Q) regimes. The alkali metals are a typical 
class of systems considered as a workbench to test the 
different theoretical approaches developed so far for the 
dynamics of the liquid state. These systems are known 
to support well-defined oscillatory modes for the density 
fluctuations (acoustic-like modes) well outside the strict 
hydrodynamic region, down to wavelengths of a few inter- 
particle distances. 

From the experimental point of view, it is worth men- 
tioning the pioneering inelastic neutron scattering (INS) 
experiment by Copley and Rowe in liquid rubidium, 
while more recently many experimental efforts have been 
performed aiming to realize more and more accurate ex- 
periments: INS investigations have been devoted to liq- 
uid cesium |^ , sodium |^ , lithium j|] , potassium ||^ and 
again rubidium |^ . Simultaneously many theoretical and 
numerical studies have also been reported on the same 
systems [7-14 . Among others, one significant advantage 
of the numerical techniques with respect to INS is the 
possibility to overcome the experimental restriction, i. e. 
the limited (Q — E) accessible region and the simulta- 
neous presence of coherent and incoherent contributions 
to the INS signal. Moreover, numerical techniques al- 
low to access a wider set of correlation functions while 
the inelastic scattering experiments basically probe the 
density-density correlation function. 



More recently, the advent of Inelastic X-ray Scatter- 
ing (IXS) technique p5|jl^ ], opened new possibilities for 
the experimental determination of the dynamic structure 
factor in the nanometer wavevector region. In fact, the 
use of X-rays often allows to extend the accessible ex- 
changed energy region to neutrons (an important merit 
especially in these systems characterized by a high value 
of the sound velocity), and to overcome the presence of 
the incoherent scattering when one is also interested in 
the collective motion. On the other hand, the resolution 
demand is particularly severe in this case, and the scat- 
tered intensity rapidly decreases on increasing the atomic 
number of the sample. This complementarity of IXS and 
INS allowed to improve the comprehension of the micro- 
scopic dynamics in alkali metals through a number of 
experiments performed in the last decade p^[2^]. 

A basic idea suggested by Molecular Dynamics (MD) 
simulations ||2^ , and recently demonstrated by IXS in the 
case of lithium |^,^ , aluminum , and sodium , 
is the existence of two different time-scales which drive 
the viscous decay of the density fluctuations. These two 
time-scales are believed to reflect two different phenom- 
ena. The slower process (whose characteristic time spans 
from few picosecond to seconds depending on the tem- 
perature) is responsible for the highly correlated atomic 
motion. Its current understanding is provided by the 
Mode Coupling Theory (MCT) approach, which gives a 
satisfactory description of a large number of experimen- 
tal and numerical results obtained in glass forming sys- 
tems This process is usually referred to as the a 
(structural) relaxation process. The faster process (char- 
acteristic time in the sub-picosecond time range) is less 
understood, and, in the normal liquid phase, is usually 
described by un- correlated hmary collisions between par- 
ticles 1^. 

Recently, theoretical (hard spheres) [E8| and numeri- 
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cal (Lennard- Jones) studies on the high frequency 
dynamics of model monatomic glasses have shown that, 
even in a harmonic glass, it exists a fast relaxation pro- 
cess. This process shows the typical phenomenology of 
relaxations, as the sound velocity increases with wavevec- 
tor Q, in spite of the fact that all the diffusive degrees 
of freedom are frozen and the anharmonicity is absent. 
This findings point towards the association of this relax- 
ation process to the topological disorder i. e. the 
dephasing of the different oscillatory components present 
in the dynamics of the density fluctuations at a given Q 
value gives rise to the decay of the correlation functions. 
If this interpretation is correct, this fast relaxation pro- 
cess should be a general property of disordered systems 
and, in particular, it must also exists in the liquid phase. 

In this Paper, we investigate whether in a simple 
monatomic liquid it is possible to associate the fast re- 
laxation process to the topological disorder and to im- 
prove our understanding of this phenomenon beyond the 
present belief based on (almost) un-correlated binary col- 
lisions. This study compares the dynamical properties of 
a monatomic system (lithium) in its liquid, supercooled, 
glassy and crystalline phases obtained by Molecular Dy- 
namics (MD) simulations. The analysis of quantities such 
as the speed of sound, the acoustic attenuation, and the 
related relaxation times indicates the existence of the 
two- relaxation processes in all disordered phases. The 
slow process, whose relaxation time largely increases ap- 
proaching the structural arrest (glass transition), effec- 
tively behaves as the a relaxation. Most importantly, we 
find that the fast process is always present and its fea- 
tures do not change between the liquid and glass phases. 
This implies the importance of the topological disorder 
in this fast dynamics even in the liquid phase. Moreover, 
it also indicates that even at the short timescale consid- 
ered the atomic motion is highly correlated and cannot 
be described by uncorrelated binary collisions. 



II. MOLECULAR DYNAMICS SIMULATION 

We numerically investigated the molecular dynamics 
of 2000 lithium atoms interacting by the Price- Tosi pseu- 
dopotential. This interaction potential has been shown 
to accurately reproduce both the structural and the dy- 
namical properties [||. Lithium has been chosen due to 
the availability of high resolution experimental data on 
the dynamic structure factor S{Q, lu) in the liquid phase, 
which show the existence of the two relaxations processes 

The molecular dynamics simulations have been per- 
formed in the microcanonical ensamble, utilizing a Ver- 
let algorithm with an integration time dt = 0.5 fs. The 
system has been equilibrated at several different temper- 
atures, in a square box with periodic boundary condi- 
tions. To make closer contact with the IXS experiments. 



the numerical densities, i.e. the box dimensions, have 
been adjusted at each temperature to the experimental 
values , as far as the mass of each atom, chosen as the 
isotopic mixture corresponding to the natural abundance 
(to — 1.15221 X 10~^^ kg). The energy conservation, all 
over the explored temperature range, was better than few 
parts in 10"''. 

The melting temperature of the selected inter-atomic 
interaction potential is T„j = 450 K, in good agreement 
with the experimental melting point. Below this tem- 
perature, the "MD liquid" can be supercooled avoiding 
the crystallization, and it attains the structural arrest at 
Tc ~ 260 K, a value extrapolated from a power law fit 
of the calculated diffusion coefhcient. This value of Tc 
is consistent with that estimated from the depth of the 
effective interaction potential well, e. In fact, using the 
relation Tc « 0.47e, valid for the Lennard Jones potential 
we derive for the Price pseudopotential Tc = 270 K. 
Below Tc it is no longer possible to equilibrate the sys- 
tem, but we can obtain the simulated "lithium glass" by 
rapid quenching down to few Kelvin. Further heating of 
the system up to (but below that) Tc allows to investigate 
the effect of the anharmonicity in the glassy configura- 
tion. We studied the system in the range from 5 K to 
1100 K, a region much wider than that accessible in the 
real system. Indeed, the low and high T limits cannot be 
reached experimentally because of the unavoidable crys- 
tallization below Tm or of the high chemical reactivity 
of the hot liquid. At any of the selected temperatures 
(above Tg), we equilibrated the system for few hundreds 
picoseconds and then we stored the configurations, ri{t), 
in order to calculate the dynamic structure factor, de- 
fined as: 



SiQ,u;) = (l/iV)FT^(e"'*5-""'(*)e''^"'j(°)). 



(1) 



For the sake of simplicity, we turned around the cal- 
culation of the correlation by making use of the Wicner- 
Kintchine theorem, that simplify the computation of the 
dynamic structure factor as: 



S{Q,u;) = - 



ftj: 



iQ-Ri(t; 



(2) 



The numerical resolution 27r/At, being At the acqui- 
sition length, was set to 1.3 ps~^, while the 'free spec- 
tral range' (i.e. the accessible frequency window), 2Tr /St 
where St is the time separation of the stored configura- 
tions, was set to 120 ps"^. Each S{Q,uj) at a fixed Q 
value was obtained as an average over different statisti- 
cally independent runs (up to 52 for each temperature) 
and Q orientations in the simulation box. 
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III. RESULTS 

In Fig. |l|, we report a comparison of the simulated 
dynamic structure factor at Q = 7 nm^^ with the cor- 
responding experimental result for those temperatures 
where data are available (T=475 K and T=600 K). 
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Dynamic structure factors of lithium at 
^ and T = 600 K (upper panel) and T = 475 



FIG. 1. 
Q = 7nm 

K (lower panel) are reported together with the experimental 
data (points with error bars) |p2|. The resolutions (FWHM) 
are: SEexp ~ 4.5 ps~^ for the IXS data, SEsim = 1-3 ps~^ 
for the simulations. In the experimental data, the detailed 
balance factor has been removed to be consistent with the 
classical character of the simulation. 

The detailed balance factor has been removed from the 
experimental data through the relation 



Ic{Q,u;) = 



hio/KT 



-huj/KT 



(3) 



where the suffix c and q indicates the classical and quail- 
tum quantities, respectively. 

The good agreement between the simulated and the 
experimental spectra testifies that the adopted potential 
represents well the lithium system Q 

In Fig. ^ we report the temperature dependence of the 
dynamical structure factor for the same selected value of 
the momentum transfer, Q = 7.0 nm^^. 

The temperature dependence of the spectra reported 
in Fig. ^ shows two apparently opposite behaviors; on one 
side the central line -or Mountain peak in the terminol- 
ogy of the glass transition phenomenology- gets narrower 



with decreasing temperature and its width -related to the 
inverse of the relaxation time r- becomes negligible at 
the glass-transition, thus indicating the structural arrest 
(i. e. a divergent relaxation time). On the contrary the 
Brillouin component shows a much reduced temperature 
dependence, and the excitation maintains a significant 
linewidth even at very low temperature and no remark- 
able behavior is observed around Tr- 
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FIG. 2. Dynamic structure factor of lithium at 

Q — 7nm~^ at different temperatures and phases: liquid, 
supercooled and glass (full line), crystal (oriented along (1 
0), dashed lines). The energy resolutions (FWHM) is 
SEsim — 1.3 ps~^ (in the crystal case it has been enhanced by 
a factor three due to the narrower lineshape) and in all cases 
it is negligible with respect to the intrinsic spectral width. 



*The experimental resolutions of the reported spectra are 
slightly different for the experiment and simulation, in par- 
ticular: Sujexp ~ 4.5 ps~^ and SQexp = 0.35 nm~^ FWHM, 
while 5LJs^m = 1.3 ps"^ and SQsim = 0.2 nm~^ FWHM. The 
energy broadening of the spectra for the considered Q- value 
is, however, always larger than any resolution effects. 



It is worth to remember that the width of the Brillouin 
peaks is a non-monotonic function of T, they attain a 
maximum in the temperature region where f2r « 1, and 
become negligible when VIt » 1 (here Q. is the frequency 
of the excitation, i. .e. approximately the Brillouin peak 
position). The origin of the non vanishing width of the 
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Brillouin peak has to be searched in a phenomenon differ- 
ent from the relexation process responsible for the struc- 
tural arrest (a-process) . 

To investigate the effect of a possible origin (the an- 
harmonicity) of the observed Brillouin linewidth, we re- 
port the spectra obtained by MD at 5 and 300 K in the 
lithium crystal, produced by thermalizing at finite tem- 
perature the configuration obtained starting from an or- 
dered (BCC) lattice. 

As clearly shown in the figure, and as well known 
from the experiments, the phonon linewidth in the crys- 
tal is completely negligible compared to the disordered 
phases, indicating that the residual linewidth cannot be 
attributed to anharmonicity, which, in turn, is expected 
to be of the same order of magnitude in the ordered and 
disordered phases. 
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FIG. 3. Dynamic structure factor of a lithium simulated 
glass at T=5 K. 

Having excluded the anharmonic origin of the Brillouin 
linewidth, a possibility lies in the presence of a further 



relaxation process, that in the system parallels the a- 
process. In order to investigate on such a possibility, an 
appropriate way can be to study the dynamic structure 
factor in the glassy phase, where any relaxation effect due 
to the a process is frozen. The dynamic structure factor 
of a lithium glass, obtained by "instantaneous" quench 
from above the melting point, is reported in Fig.^ 

As a demonstration of the presence of such an addi- 
tional relaxation process, we report the observation of a 
positive dispersion of the sound velocity below the melt- 
ing temperature. This effect, also expected to take place 
at frequencies such that ujt « 1, is shown in Fig. ^, where 
we report the quantity ci{Q) ~ '-^i{Q)/Q with uji the max- 
imum of the current correlation spectra 



J{Q,u;)=u;yQ^SiQ,Lu). 



(4) 



The existence of a positive dispersion of sound, repro- 
duced by the simulations of the glass and cold liquid 
phases (it is missing in the crystal), supports the hy- 
pothesis of the existence of a second relaxation process. 
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FIG. 4. Speed of sound in different thermodynamic 
phases. Dots:glass/supercooled, downtriangle; crystal along 
(1 0), uptriangle: crystal along (1 1 1). Full and open sym- 
bols are relative to T = 5 K and T = 300 K, respectively. 
Arrows indicates acoustic data at T = 70 K. Lines are guide- 
line for the eyes only. 

The qualitative analysis of Fig. ^ per se indicates the 
presence of at least two relaxation processes governing 
the high frequency collective dynamics of the consid- 
ered system: one, strongly temperature dependent, which 
shows up mostly in the width of the central line, and 
the other, strongly temperature independent, which gives 
rise to the width of the inelastic feature. These simu- 
lations confirm, therefore, in a much wider temperature 
range the conclusions provided by the analysis of the IXS 
data on lithium [ p2[ , and extend into the supercooled 
and glassy phase regions the theoretical predictions of 
a dynamics ruled by two different time-scales. Partic- 
ularly interesting is the persistence of the second (fast) 
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relaxation process, which is essentially identical in the 
liquid and in the glass. This result is somewhat in con- 
tradiction with the current belief: in fact, in the liquid 
the fast dynamical processes has been so far associated 
to uncorrelated binary collisions, while in the solid the 
atomic motion is dominated by long range correlations. 
A quantitative description of this relaxation process may 
contribute to shed light on this point. 
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FIG. 5. Seletecd fitting examples at different Q—T values. 
Open dots: molecular dynamics, fuUlines: fits. ^(Q, 0) values 
at T = 5 K are 2.1 x 10"^ and 4.9 x 10"^, respectively. 

We now quantitatively represent the dynamic structure 
factor to asses the evolution of the relaxation processes 
from the liquid down to the supercooled and glassy state. 
To this purpose, we analyzed our data adopting a model 
based on the memory function approach, introduced to 
represent the MD data on alkali metals |3^j3^, and re- 
cently utilized to interpret the experimental data of liquid 
lithium ||2^ , aluminum and sodium . Within this 
framework, the evolution of the density autocorrelation 
function is described by a memory function trough the 
Langevin equation: 

(g, t) + ujl{Q)(f>{Q, t) + Jo Af(Q, t-t')'(l) (g, t')dt' = 



where 



(5) 



By Fourier transformation, the dynamic structure fac- 
tor is given by: 



S{Q,oj) ^ TT-^uol{Q)M'{Q,uj) 

S{Q) ^ [t^2 _ ^2 _^ ujM"{Q, u)f + [ujM'iQ, uj)y 



(6) 



The complex memory M(Q, t) contains all the inter- 
action details and decays over the characteristic time- 
scales of the system. The dynamics of a simple liquid 
can be satisfactorily described by the two time exponen- 



tials ansatz |33 



AfL(g,i) = A2(g)e~*/^=(e)-f A2(g)e- 



-t/^M(Q) 



(7) 



where the indexes a and /i indicate the usual a- (or 
structural-) and the fast- (or microscopic) process, re- 
spectively. We used Eq. (||) as a model function. Its con- 
volution with the simulation resolution R{uj) has been 
utilized as a fitting function to the MD spectra: 



F{Q,uj) = J S{Q,Lu')R{uj-Lu')dLj' 



(8) 



The results of such procedure are shown in Fig. 
where we report the calculated dynamic structure fac- 
tors together with their best fitted line-shapes at selected 
momenta and temperatures. 
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FIG. 6. Q-Dependence of the relaxation times at some 
selected temperatures. The slight increase of the a process at 
low Q is an artifact due to finite resolution effects (see text). 

Among the fitting parameters, the relaxation times of 
the two processes are reported in Fig. ^ as function of 
the exchanged wavevector and at selected temperatures. 
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At all considered temperatures, the a-process shows an 
almost flat Q dependence, with an apparent increase be- 
low Q ~ 4 nm^^, an artifact that occurs when the re- 
laxation time is comparable or higher than the inverse 
resolution, already observed in the analysis of the IXS 
spectra |22 . The microscopic process shows, instead, an 
effective slightly decreasing Q trend. The value of the 
structural relaxation time at large Q shows a sharp in- 
crease in the supercooled region and upon further cool- 
ing. On the contrary, the temperature behavior of the 
fast microscopic, /x, process is very different. In fact, it 
gives a basically T-independent relaxation time in any of 
the considered phases. This figure provides a quantita- 
tive account of the qualitative analysis that was derived 
from the results reported in Fig. 1. 
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FIG. 7. T-Dependence of the relaxation times (circles, left 
axis). In the same plot we also report the behavior of the 
diffusion coefficient (stars, right axis). 

To be more precise in the assessment of the tempera- 
ture dependence of the relaxation times, once observed 
that the Q— dependence is similar for each tempera- 
ture, we defined a Q-independent temperature coefficient, 
7(T), through the relation t(Q,T) = "f{T)T{Q,To). As 
arbitrary reference temperature we selected To = 1100 K, 
and we determined 7(T) as the scaling factor that min- 
imizes the mean squared differences between each set of 
data and the set at To. 

The result of such a procedure, for both the a and fi 
processes, is reported in Fig. |^. As clearly seen, the re- 
laxation time of the structural process increases about six 
times between Tq and the lowest available temperature 
(T — 275 K). In the same plot we also report the quan- 
tity T/D, being D the diffusion coefficient determined 
in the simulations. As expected, this quantity follows a 
behavior similar to that of the relaxation time of the a- 
process, thus confirming the validity of its identification 
with the structural relaxation process. On the contrary, 
the microscopic relaxation time exhibits a small varia- 



tion in the whole temperature range with a total change 
of only 50% of the value at To, and without any notewor- 
thy behavior around Tc- This demonstrates that the fast 
relaxation process is the same from the hot liquid down 
to the glass. 



IV. DISCUSSION AND CONCLUSIONS 

In the present paper, we have shown how the dynamics 
of a simple monatomic system can be described in terms 
of two viscous relaxation processes in a range of tempera- 
ture covering all disordered phases: the hot liquid, the su- 
percooled liquid and the glass. Few comments concerning 
the physical origin of these two processes are in order. By 
its strong temperature dependence and by the close rela- 
tion with the temperature dependence of the diffusion co- 
efficient, the slower relaxation process is unambiguously 
identified as the a process, i. e. the process that drives 
the glass transition in those systems capable to sustain 
strong under-cooling. More interesting are the indica- 
tions obtained from the analysis of the faster process. 
The comparison between glass and crystal, the equiva- 
lence of the fast relaxation process in the liquid and in 
the glass, and the consistency between the simulation and 
experimental results in the liquid, are all indications that 
the origin of this process must be searched into something 
that is common between liquid and glass, but not existent 
in the crystal. On the considered fast time scale, which 
corresponds to values of fractions of ps, one can con- 
sider all diffusion processes in the liquid frozen, and on 
this snap-shot time-scale one can consider that the liquid 
and the glass are not only topologically, but also dynam- 
ically equivalent. Along this line, one is then driven to 
conclude that -similarly to what has been found in other 
model glasses also in the liquid this fast relaxation 
process is intimately related to the topologically disorder, 
while anharmonicity and/or dynamical effects play only 
a secondary role. The way in which the disorder can pro- 
duce the phenomenology of a "relaxation process" -also 
in the harmonic limit- is as follows. On a general ground, 
a relaxation process can be pictured as the macroscopic 
manifestation of microscopic phenomena associated with 
the existence of channels by which the energy stored in 
a given "mode" relaxes towards other degrees of free- 
dom. The spectrum of the density fiuctuations, S'((5,w), 
through the fluctuation-dissipation relation, reflects the 
time evolution of the energy initially stored {t=to) in a 
plane wave of wavelength 2it /Q. As the plane wave is not 
an eigenstate of the disordered system, aXt>to there will 
be a transfer of amplitude from this plane wave towards 
other plane waves with different Q values. This process is 
controlled by the difference between the considered plane 
wave and the true normal modes of the topologically dis- 
ordered glassy structure. This energy flow takes place on 
a specific time scale r, and gives rise to the observed re- 
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laxation process phenomenology. The depicted scenario 
is in contrast with a description of the fast relaxation pro- 
cess in terms of (more or less correlated) binary collisions. 
Indeed, according to the usual description, the collision 
between two atoms -at a given position in space f* and 
at a given time t*- are not correlated with the dynamics 
of the atoms far away from f* , and with the dynamics of 
the colliding atoms at time larger than t* + t^. In this 
view acts a the correlation time of the local dynamics. 
On the contrary, according to the present description, the 
dynamic of the whole system is strongly correlated -and 
described in term of normal modes of vibration- up to a 
time of the order of Tq, , and indicates the decorrelation 
time for the specific plane waves, not for the dynamics. 

The common origin of the fast process in the liquid and 
in the glass, and its association to the disorder, poses 
against its interpretation in terms of kinetic processes 
such as "binary collisions" . In this sense, the idea that 
all the binary processes are contained in the short time 
dynamics (THz region), while the correlated motion is 
fully described by the a process has to be revised. 
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